DEPHASING DUE TO NONSTATIONARY 1/f NOISE 
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Motivated by recent experiments with Josephson qubits we propose a new phenomenological 
model for noise due to collective excitations of interacting defects in the qubit's environ- 
ment. At very low temperatures the effective dynamics of these collective modes are very slow 
leading to pronounced non-Gaussian features and nonstationarity of the noise. We analyze the 
influence of this noise on the dynamics of a qubit in various regimes and at different operation 
points. Remarkable predictions are absolute time dependences of a critical coupling and of 
dephasing in the strong coupling regime. 

Despite a large number of investigations, the detailed explanation of the origin of 1/f noise 
stays an open problem in condensed matter physics. The diversity of systems that exhibit 1/f 
fluctuations suggests that its physical origin is not universal^ Still, because of its ubiquity, the 
study of underlying mechanisms, relevant at least in some cases, is interesting. 

In the context of quantum-information processing the field has attracted new attention. In 
particular in solid state devices, e.g. Josephson qubits, low frequency noise, typically with a 1// 
power spectrum, is considered as the most important limitation for the preservation of the phase 
coherence. The microscopic origin of this noise depends thereby on the considered system. In 
the case of charge noise, for instance, it is believed to be due to the activity of random traps for 
charges in the dielectric substrate and in the oxide barrier of the Josephson junctiodP Similarly, 
neutral defects in the junction barriers can cause critical current fluctuations and low- frequency 
flux noise may take its origin from trapped vortices and magnetic impurities. In practice, 1/f 
noise appears difficult to suppress and, since dephasing is dominated by low-frequency noise, it 
is particularly destructive. 

Decoherence of a qubit due to 1/f noise has already been studied in literature using various 
models for the noise: One possibility is to assume Gaussian-distributed fluctuations, e.g. of 
a large collection of microscopic modes, which in some cases may be modeled by a harmonic 
oscillator bathP Alternatively, one can consider a microscopic model with independent bistable 
fluctuators, each one characterized by a switching rate 7^ and an effective bias voltage V{ induced 
at the qubitP^If the switching rates are distributed according to ^(7) ~ I/7 the resulting power 
spectrum exhibits a 1/f divergence at low frequencies. If only a few fluctuators (microscopic 
modes) dominate, the effect of the noise differs. 

In this work, we present a new model of 1/f noise and discuss its effects on the phase 
coherence of a qubit. In contrast to the mentioned proposals we consider situations where 
the mutual interaction of the localized defects in the background cannot be neglected, e.g. at 
very low temperatures. Since neutral and charged defects behave as elastic and electric dipoles 
respectively^ the strength of the interaction is long-range in both cases decaying as 1 /r 3 as 
a function of the distance r. Although weak, this interaction may play a crucial role at ver 
low temperatures, leading to formation of highly cooperative clusters of interacting defect; 
The resulting slow dynamics can be described as a perpetual search of the absolute free-energy 
minima: the system spends most of the time in local free energy minima corresponding to 
metastable configurations, occasionally flipping between these states. The exponential relation 
between the typical trapping time r and the energy barrier to overcome E, t(E) ~ TQe E ^ T , 
implies that the distribution of the trapping times, P(r) is generally broad, decaying very 
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slowly for large times: P(t) — r _1_M . Here \x is a positive dimensionless number that can be 
written as the ratio of temperature T and a typical height i?o of the energy barriers, fx ~ T/Eq. 

The fluctuating extra bias induced at the qubit is described by an asymmetric telegraph 
noise, where the times between two successive plateaux are distributed according to an algebraic 
distribution 

The heights and lengths of the plateaux are distributed according to narrow distributions V(v), 
which is assumed to have zero mean, and p(r|) respectively. The coupling strength to the qubit 
is characterized by the dimensionless parameter g 2 = ((utj) 2 ). For [i > 2, when r = drrP(r) 
and t 2 are finite our model is equivalent to usual Poissonian telegraph noise as proposed u¥^l 
However, at ultra-low temperatures, i.e. fi < 2, important differences emerge." We thereby have 
to distinguish two additional classes of noise: (i) for 1 < \i < 2, r is still finite but r 2 diverges. 
Consequently, the rate of occurring plateaux fluctuates strongly around the average r _1 . In 
particular, due to the diverging r 2 the statistics of the noise are only approximately Gaussian, 
even in the weak coupling regime, (ii) for < \i < 1, when both r and r 2 diverge, it is even 
impossible to define an average rate of occurring plateaux, i.e. the noise is scale invariant. In 
this case, the two point noise correlator decays as 



x(t + t)x(t) - x(t + 1>) x(t) ~ x(t) {t/t'y-», (2) 

as opposed to the exponential decay in the case \x > 2. For < fj, < 1, the power spectrum thus 
exhibits a l//'* divergence at low frequencies. Furthermore, the t/t'-scaling in © implies an 
intrinsic nonstationarity of this noise. 

In this work we focus on the influence of such a noise on the phase coherence of a qubit. 
The full Hamiltonian of a dissipative quantum two-level system can be written in the form: 

H = --e a z - - (smrja x + cost/ a z ) X + H cnv (X) , (3) 

where e and r\ are control parameters which we assume to be time-independent and X is the 
fluctuating extra bias induced by the environment described by 7i erw . For definiteness we con- 
sider in the following two special working points for the qubit: (i) the optimal point (rj = ir/2), 
where the linear longitudinal coupling to the noise is tuned to zero and consequently decoher- 
ence is reduced considerably, ^ and (ii) the case of longitudinal coupling (rj = 0). Below we will 
show that transverse linear coupling is - regardless of the considered noise model - equivalent 
to quadratic longitudinal coupling and a remaining transverse noise component with suppressed 
low-frequency tail. 

Optimal point (r/ = 7r/2): In general, transverse noise leads to relaxation, but also contributes 
to pure dephasing in higher orders. The effect of low-frequency transverse noise (uj <C e) can 
be treated in the adiabatic approximation: the state of the qubit follows the effective magnetic 
field and © can be diagonalized to -o z \/e 2 + X 2 /2 ps -a z [e + X 2 /(2e)] /2. Therefore, low- 
frequency transverse noise contributes to pure dephasing as quadratic longitudinal noise. On 
the other hand, higher frequencies, uj ~ e, mainly contribute to relaxation and the adiabatic ap- 
proximation breaks down. Situations where both low and high frequencies are present therefore 
need further analysis. 

The unitary transformation U(t) = exp{iQ(t)a y /2} with O(i) = arctan(X(i)/e) transforms 
the Hamiltonian (j2J) to a spin frame that follows the effective field ez + X(t)x: 

H = UWUHt)+iU(t )±uHt) = -\^e 2 + X 2 (t) r z + l -Q{t) r y , (4) 



"Our predictions also differ from those for an ensemble of independent Poissonian fluctuator jffi^ A detailed 
comparison is not straightforward and will be addressed elsewhere. 



where the Pauli matrices Tj remind that the spin is now measured in the time-dependent basis, 
defined by r z (t) = U(t)a z U'(t). Note that for small amplitude of the noise |-X"(i)| <C e the new 
basis is almost identical to the original basis. 

As we can see from the expression of the transformed Hamiltonian (J3J) linear transverse 
noise is indeed equivalent to quadratic longitudinal noise - as predicted in the adiabatic ap- 
proximation. But due to the explicit time dependence of the transformation U(t) an addi- 
tional transverse component proportional to the time derivative of the original noise X appears, 
6(t) = (X(t)/e)/{l + (X(t)/e) 2 ) ps X(t)/e. One might think that the transformed Hamiltonian 
is therefore even more complicated than the original one since we are now dealing with both 
longitudinal and transverse fluctuations. But a closer look at the expression reveals that the 
low-frequency transverse noise is suppressed, Sq(uj) = 2(Q(t)Q(t')) (JJ ps (^) 2 Sx(w). Hence, in 
the experimentally relevant situation where the spectrum Sx is less singular than 1/uj 2 at low 
frequencies, Sq is a regular function at low frequencies with Sq(u> = 0) = 0. As a consequence 
the transformed transverse noise does not contribute considerably to dephasing but only to re- 
laxation. Dephasing is dominated by the longitudinal part of the noise proportional to X 2 that 
still contains all low-frequency contributions. We can therefore add up the two contributions in 
PJ independently and extract simple expressions for dephasing and relaxation rates: 

Tf 1 = ±S e (u; = e)^\sx(u = e) (5) 

T 2 l = \T^+T- 1 (6) 

where T\ and the first contribution to Ti are due to the second term in (j3J) and can be calculated 
simply using the golden rule. The pure dephasing time T v is defined as the characteristic decay 
time due to quadratic longitudinal noise X 2 (t)/(2e) in In the case of a non-singular noise 
spectrum it is given by T" 1 = S x z(u = 0)/(4e 2 ). For Gaussian 1// noise it has been calculated 
m RefP In this case the dephasing law differs from a simple exponential and one should add two 
decay laws rather than the rates in ©. For the model of nonstationary 1// noise considered in 
this work the evaluation of T v can be reduced straightforwardly to the case of linear longitudinal 
coupling: X(t) and X 2 (t) /(2e) both describe an asymmetric telegraph noise, the major difference 
being a change in the amplitude and thus in the coupling strength. 6 The coupling strength at 
the optimal point is thereby considerably reduced in comparison to linear longitudinal coupling. 

Longitudinal coupling (rj = 0): We now focus on the case of pure dephasing, r\ = in ©• 
Therefore, we follow the time evolution of the off-diagonal entries of the qubit's density matrix, 
(cr+(i)) = D(t p ,t)e ltt (cr+(0)) , where the dephasing factor is defined by D(t p ,t) = e**^'*) with 
<fc(tp,t) = ftp +t X(t')dt' . Because of the intrinsic nonstationarity of the noise © we expect 
to find an explicit dependence of dephasing on the preparation time t p . In the following we 
will therefore investigate this explicit dependence of D(t p ,t) on t p . This can be achieved using 
CTRW techniques and taking special care of initial conditions at tjr^ 

As already discussed previously we have to distinguish three different cases depending on 
the value of //: When fi > 2 the noise is equivalent to Poissonian telegraph nois and therefore 
stationary (for t p > To). Consequently, dephasing is insensitive to the absolute time t p . 

For 1 < /j, < 2 it turns out that nonstationarity only plays a role in the strong coupling 
regime g > 1. For weak coupling g < 1 the noise is qualitatively equivalent to the Poissonian 
case fi > 2, i.e. only subdominant corrections appear^! For g > 1 the dephasing law differs 
from a simple exponential and the dephasing time, defined as D(t p ,T^,) = e _1 , scales as ~ 
To[l/e + (to/(to + ip))^ 1 ] -1 /^ -1 ). The dependence of on t p thus only matters for values of 
/i close to 1 and disappears as fi increases to higher values (see Figure P). 

^In addition X 2 (t)/(2e) has a finite mean. A thorough discussion of this case will be presented elsewhere. 



Figure 1: Dephasing time as a function of the coupling strength g for /i < 1 (left) and /i > 1 (right). 



As expected nonstationarity turns out to be most prominent in the case < fj, < 1. Again, 
we have to distinguish two regimes of weak and strong coupling separated by a critical coupling 
constant g c . But here g c is not of order one, but depends explicitly on the age of the noise: 
g c {tp) — A(/i)t p ^ 2 <C 1- Therefore, the range of the strong coupling regime increases with the 
age of the system implying that any qubit subject to a noise with < fi < 1 eventually ends up 
in the strong coupling regime. 

For weak coupling g < g c (t p ) dephasing is essentially insensitive to the age of the noise. Away 
from a (short) initial aging regime D(t p ,t) decreases as exp(— (t/r^)^), where the dephasing 
time scales as ta, ~ To<?~ 2 /^ (in the marginal case /i = 1 we get logarithmic corrections ~ 
Tog~ 2 \lng\). Indeed, for t < r^, D(t p ,t) is accurately given within a Gaussian approximation 
of the phase: D(t p ,t) ~ exp(— ^(t p , t) 2 ). Only for t > the Gaussian approximation breaks 
down and the exponential decay goes over to a much slower algebraic one, D(t p ,t) oc 

In the strong coupling regime g > g c (t p ) the dephasing depends strongly on the age of the 
noise. The Gaussian approximation of the phase breaks down even at short times and the decay 
is not exponential any more. The dephasing time becomes proportional to t p , ~ C(iJ,)t p , and 
saturates as a function of g. We emphasize that for \i — > 1, C(/i) takes smaller and smaller values 
implying that the dephasing time can be orders of magnitude shorter than the preparation time. 

In summary, we have proposed a new model of noise describing phenomenologically 

a slow collective environment. The present model shows a cross-over from nonstationary 
noise to a memoryless Poissonian telegraph noise as a function of temperature. We analyzed 
the decay of coherence of a qubit subject to this noise at different operation points. For low 
temperatures and strong coupling the dephasing time was found to depend explicitly on the age 
of the noise. 
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